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Abstract 



This article is the continued version of the analytical blowup solutions for 2-dimensional 

CN ■ 

^ ■ Euler-Poisson equations [12] , With extension of the blowup solutions with radial symmetry 

■ for the isothermal Euler-Poisson equations in , other special blowup solutions in with 



non-radial symmetry are constructed by the separation method. We notice that the results 
are the first evolutionary solutions with non-radial symmetry for the system. 
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X ' 1 Introduction 

The evolution of a self-gravitating fluid (gaseous stars) can be formulated by the isentropic Euler- 
Poisson equations of the following form: 

Pt+V ■ (pu) =0, 

{pu)^+y ■ {pu (g) u) + VP =-pV$, (1) 

A$(t,f) = a{N)p, 

where a (TV) is a constant related to the unit ball in i?^: a(l) — 2; a (2) — 2tt and For iV > 3, 

a{N) = N{N - 2)V{N) ^ N{N ~ 2) ^^^^^ ^ , (2) 
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where V{N) is the volume of the unit ball in and T is the Gamma function. And as usual, 
p — p{t, x) and u — u{t, x) — [ui, U2, un) G are the density and the velocity respectively. 

In the above system, the self-gravitational potential field $ = ^{t, x) is determined by the 
density p through the Poisson equation. 

The equation ([IJ3 is the Poisson equation through which the gravitational potential is deter- 
mined by the density distribution of the density itself. Thus, we call the system ([T|) the Euler- 
Poisson equations. The equations can be viewed as a prefect gas model. The function P = Pip) 
is the pressure. The 7-law can be applied on the pressure P{p), i.e. 

P(p) = i^p^:=^, (3) 

7 

which is a commonly the hypothesis. The constant 7 = cp/cy > 1, where cp, Cy are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in ([3]). In particular, the fluid is called isothermal if 
7 = 1. It can be used for constructing models with non-degenerate isothermal cores, which have a 
role in connection with the so-called Schonberg-Chandrasekhar limit [5j . 

The system can be rewritten as 



Pf 4- V • -t- V p ■ u =0, 
P i'-ut + Eti u,§^) + £-Pip) =-Pi&Hp), for ^ = 1, 2, ...iV, (4) 



at ' ^k=l "-l-^dxk I ^ dx,^ yi-'! ^ I'dx, 

A$(t,a;) = a{N)p. 

For iV = 3, (|4]) is a classical (non-relativistic) description of a galaxy, in astrophysics. See [2], [3] 
and [8] for a detail about the system. 

For the local existence results about the system were shown in [9], [1] and [5]. In particular, 
the radially symmetric solutions can be expressed by 

p{t, x) = p{t, r) and u{t, x) = -V{t, r) := -V, (5) 

r r 

1/2 

where the radial diameter r :— ^X^ili ^^l) • Historically in astrophysics, Goldreich and We- 
ber constructed the analytical blowup (collapsing) solutions of the 3-dimensional Euler-Poisson 
equations for 7 = 4/3 for the non-rotating gas spheres [6]. After that, Makino [10] obtained the 
rigorously mathematical proof of the existence of such kind of blowup solutions. Besides, Deng, 
Xiang and Yang extended the above blowup solutions in {N > 3) [4]. Then, Yuen obtained 
the blowup solutions in R^ with 7 = 1 by a new transformation [12j . The family of the analytical 
solutions are rewritten as 
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For TV > 3 and 7 = (2iV - 2)/iV, in g] 

1 



0, for a(<)Z^ < r. "^^^ 



i(<)= 



-A 



a(0) = 7^ 0, d(0)== as, 



(6) 



(2j^^-^2)g ^(^)'^^^'^'' - 2/(0) - « > 0, y(0) = 0, 

where ^jl = [7V(7V - 2)\]/{2N - 2)^ and the finite is the first zero of y{z); 
For TV = 2 and 7 = 1, in [12] 

p(i,r) = ^e«('-M*)), 1/(^,0 

a(t)=— — , a(0) = a, > 0, a(0)= 02; 
a(t) 

a(2) 



(7) 



2/(2:)+-2/(^)H — j7-e 

Z K 



A*, y(0) = a, m = 0, 



where K > 0, fi = 2X/K with a sufficiently small A and a are constants. 

And for other special blowup solutions, the readers may see the details in [TT] . 

Very recently, Yuen extended the above solutions to the pressureless Navier-Stokes-Poisson 
equations with density-dependent viscosity in |13| ; 



TV- 1 

Pt + Vpr + pVr + pV ^ 0, 

r 

p{Vt+VVr) + fp[t,s)s^-Hs =[np% {li^V + Vr) + {Kp'Wr, 



TV- 1, 



'-Vr 



(8) 



However, the known solutions are all in radial symmetry. In this paper, we are able to obtain the 
similar results to the non-radially symmetric cases for the 2-dimensional Eulcr-Poisson equations 
(HI) in the following theorem: 



Theorem 1 For the 2-dimensional isothermal Euler-Poisson equations Q), there exists a family 
of solutions, 



p{t,x,y) = ^e 



-c 



-^1 \ ^(i) I \ 
i[t,x,y) = —-{x,y} 



a{t) 



a{t) — ai + a2t, 



(9) 



$(s) - e*e~-fe^ = 0, $(0) = a, $(0) = /3, 



where A, B, ai ^ 0, 02 , = e* > with that A and B are not both 0, a and (3 are constants. 

In particular, ai > and 02 < 0, the solutions (Qi blow up in the finite time T — —a^jai. 



Remark 2 The solutions 0) is a line source or sink in terms of cylindrical coordinates, when we 
take z direction to lie alone the characteristic line of the source or sink. For the physical significance 
of such kind of solutions, the interested readers may refer. P. 409-410 of 111 for details. 
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2 Separable Blowup Solutions 



Before presenting the proof of Theorem [U we prepare some lemmas first. 

Lemma 3 For the continuity equation in , there exist solutions, 

j Ax + By \ 

Pit,x,y) = ^ , u{t,x,y) = ^ix,y), 

a'^[t) a(t) 

where the scalar function f{s) > G and a{t) 7^ G C^. 
Proof. We plug the solutions (O into the continuity equation 



(10) 



pt + V ■ up + Vp ■ u 



f Ax + By 
H a{t) 



d_ 

dt 



a{tY 



-2a(t)^ ( Ax + By 



a{tf ' V aW / a{tf dt 
'Ax + By 



f Ax + By \ f Ax + By 

a{ty a{t)^ a{tf 



1 d , fAx + By\ ait) f d d 



a{t) 



a{t) 



d 



f 



ait) 



d 



f 



dx 



X 



a{t) J ' a{t) \dx 
Ax + By" 



a{t) 
Ax + By 



a{t) 



dy a{tr^ 



y 



2d{t) ^ f Ax + By 



a(t)3 ' \ a{t) 
'Ax + By 



a{t) 
0. 



/ 



ait) 



Ax 



<t)'^[ ait) 
'Ax + By 
ait) 



f 



[Ax + By 



1 ■ ( Ax + By\ iAx + By)ait) ait) \ ait) 



a(t)2 



ait) a(t)2 



/ 



By 



aity^ 



ait) 



aitf 



ait) 



(11) 
(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



The proof is completed. ■ 

On the other hand, we may use the fixed point theorem to show the local existence of the 
following ordinary differential equation: 



Lemma 4 There exists a sufficiently small xq > 0, such that the equation 

<E.(s) 

*(s) -e*e" if = 0, 
$(0) = a, $(0) = /3, 

where e* > 0, a, and (3 are constants, has a solution $ = "J'(s) G C^[0, Sq]. 
Proof. We integrate the equation (fT5|) once: 



(18) 



$(s) 



Ke*e T^d-q + (3 



(19) 
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And set 



Ke e dri + (3. 



then for any sq > 0, we get / £ C^[0, sq]. and for any $2 £ C^iO, sq], we have, 



\f{s,<i>i{s))-f{s,<i>2{s))\^K, 



■^1 (v1 _±2M_. 

(e — e ^ )dr] 



(20) 



(21) 



As is a C"'^ function of y, we can show that the function e^, is Lipschitz-continuous. And we get, 
|/(s,$i(s))-/(s,$2(s))| =0(1) / |($i(?7)- $2(77)1 rfr,<0(l)so sup |$i(s) - $2(s)| • (22) 

Jo 0<s<so 



We let 



We have T$ e C[0, so] and 



T$(s) =a+ / /(??,$(ry))d??. 
Jo 



(23) 



lT$i(s)-T$2(s)| 



-/o 



<0(l)so sup |$(s)i -$(5)2! 

0<s<so 



(24) 



By choosing sq > to be a sufficiently small number, such that O(l)so < 1, this shows that the 
mapping T : C[0, sq] C[0, sq], is a contraction with the sup-norm. By the fixed point theorem, 
there exists a unique $(s) e C[0, sq], such that T$(s) = $(s)- The proof is completed. ■ 

And we need another lemma to show the global existence of the ordinary differential equation 
(HHI. 



Lemma 5 The equation, 

$(s) -e*e~ K =0^ 
$(0) = a, $(0) = /3, 

where e* > 0, a anc? /? are constants, has a solution in (—00, +00) and lim $(5) 

S — >±C30 

Proof. We prove the case $(0) = first. By integrating we have, 

$(s) = e* / e~ K dr]>0. 
Jo 

Thus, for < s < So, 'l'(s) has a uniform lower bound 

$(s) > $(0) = a. 



(25) 



(26) 



(27) 



As we obtained he local existence in Lemma there are two possibilities: 

(1) <i>(s) only exists in some finite interval [0, sq]: (la) lim $(s) = +00; (lb)<i>(s) has an uniformly 
upper bound, i.e. <i>(s) < uq for some constant aQ. 

(2) $(s) exists in [0, +00): (2a) lim $(s) = +00; (2b)<i>(s) has an uniformly lower bound, i.e. 

s— )--t-oo 

<I>(s) < a for some constant ao- 
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We claim that possibility (1) does not exist. We need to reject (lb) first: If the statement (lb) is 
true, becomes 

ao f 

e*se~T< >e* e~~dr] = $(s). (28) 

Thus, $(5) is bounded in [0, sq\. Therefore, we can use the fixed point theorem again to obtain a 
large domain of existence, such that [0, sq + (5] for some positive number d. There is a contradiction. 
Therefore, (lb) is rejected. 

Next, we do not accept (la) because of the following reason: It is impossible that lim <J>(s) = +00, 

s—tso- 

as from p8|. ^{x) has a upper bound in [0, sq]: 

e*se"lf > $(s). (29) 



Thus, (1251) becomes, 



$(so) = $(0) + e* / $(77)^77 (30) 
Jo 

<a + e* se K d-q (31) 

= a + e*sle-^. (32) 



Since $(s) is bounded upper in [0, sq], it contracts the statement (la), such that lim <I>(s) = +00. 
So, we can exclude the possibility (1). 

We claim that the possibility (2b) does not exist. It is because 

^{s)^t* I K dr]>e* e' K drj = e*e^ K s. (33) 
Jo Jq 

Then, we have, 

$(s) > a + e*e"^s. (34) 
By letting s +00, ((34| turns out to be, 

$(s) = +00. (35) 



Since a contradiction is established, we exclude the possibility (2b). Thus, the equation ([25]) 
exists in [0, +00) and lim <&(s) = +00. Due to the solution is symmetric about s — 0, we have 

s — >+oo 

lim <i>(s) = +00. 

s — > — 00 

For the case <i>(0) = j3, these exist constants s and ag, such that $(0) — (3 m the ordinary 
differential equation: 

$(S) = as, $(5) = 0. 
The above fact is due to the transformation z = s — s to have 

<i>(z)-e*e- K ^0, ^3^^ 
$(0) = as, $(0) = 0. 
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Then we get lim ^{z) = lim $(s) = +00 with the previous result and the continuity of the 

z — *±oo s — >±oo 

solutions ([M]) to show the lemma is true. This completes the proof. ■ 

On the other hand, the following lemma handles the Poisson equation for our solutions 



Lemma 6 The solutions, 



a{ty 



+c 



with the second- order ordinary differential equation: 



4.(s) 



*(s) - e*e~~ = 0, $(0) = a, *(0) = /3, 



(38) 



(39) 



where s := {Ax + By) / a{t) and C , = e* > 0, a and (3 are constants, 

fit into the Poisson equation ([^3 in B? . 

Proof. We check that our potential function ^{t,x,y) satisfies the Poisson equation 



A^{t,x,y) - 2'Kp 

'Ax + By 



2tt 



a(t) J a{ty 



-c 



V- 

d_ 

dx 



^ Ax + By \ A ^ fAx + By\ B 



[t] J a{t)' \ a{t) J a{t) 



ZTT \ a(t) ) 



+c 



a{tY 



^, Ax^By A 
ait) a{t) 



d_ 

dy 



$( 



Ax^By B 
ait) ''^ 



2n \ a(t) ) 



+ C 



a{t)^ 



A^ + B^ 
~^(tj^ 



2^gC 



A^ + B^ 



re K 



where A and B are not both 0. 

Then, we choose s :— (Ax + By)/a{t) and the ordinary differential equation: 

$(s) - e*e'~ = 0, $(0) = a, $(0) = /?, 



(40) 
(41) 
(42) 
(43) 
(44) 



(45) 



with ^2^Q2 — £* , a and f3 are constants in Lemmas S] and [SJ Therefore, our solutions ([55]) satisfy 
the Poisson equation 

The proof is completed. ■ 

Now, we are ready to check that the solutions fit into the Euler-Poisson equations (jl]) . 
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Proof of Theorem [TJ By Lemma |3] and Lemma [51 the solutions (O satisfy and For 
the x-component of the isothermal momentum equations (|3])2 in i?^, we have 



dui 
~dt 



dui 



dui 



h U2- „ 

ox ay 



dx dx 



d_ (a{i^ \ ait)_ d_ / aft) 

[ Ax±By\ 



a{t) d_ f djt) 



K 

- P 
-0, 



d e 



-+c 



dx 

a{ty 



a{tf 

- p<i> 



Ax^By\ A 



i{t) 



Ax + By\ A 



a{t) J a{t) 



p^ 



a{t) 

I'Ax + By'X A 
V ait) J W) 



by taking d{t) — that is 



a{t) — ai + a2t. 



(46) 
(47) 

(48) 
(49) 
(50) 

(51) 



For the y-component of the isothermal momentum equations the proof is similar. As the 

readers may check it, we omit the detail here. 

We have shown that the solutions ^ satisfy the Euler-Poisson equations. In particular, oi > 
and 02 < 0, the solutions ([9]) blow up in the finite time T = —02/01. 
The proof is completed. ■ 

Remark 7 For the case of A — and B ~ 0, the corresponding solutions 0) may be deduced to 
the special solutions in ill}/ . 

Remark 8 Our solutions (0i also work for the isothermal Navier-Stokes-Poisson equations in : 



Pt+V ■ [pu) =0, 
(pu)j+V ■ {pu®u) + VKp =-pV$ + pl^u, 
A$(t,a;) = 27rp, 



(52) 



where n > is a positive constant. 



Additionally, the blowup rate about the solutions is immediately followed: 



Corollary 9 The blowup rate of the solutions (0) is, 



limp(t,0,0) {T-ty > 0(1). 



(53) 



In conclusion, due to the novel solutions obtained by the separation method, the author conjec- 
tures there exists other analytical solution in non-radial symmetry. Further works will be continued 
for seeking more particular solutions to understand the nature of the Euler-Poisson equations (j4]) . 
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